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Abstract: Photons produced in heavy ion collisions escape virtually unperturbed from 
the surrounding medium, thus representing an excellent probe of the conditions at the 
emission point. Using the gauge/gravity duality, we calculate the rate of photon production 
in an anisotropic, strongly coupled M = 4 plasma with N { <C N c quark flavors. We consider 
arbitrary orientations of the photon momentum with respect to the anisotropic direction, 
as well as arbitrary values of the anisotropy. We present results for the correlation functions 
of two electromagnetic currents and for the electric conductivity. These quantities can be 
larger or smaller than the isotropic ones, depending on the direction of propagation and 
polarization of the photons. 
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1 Introduction 

As the photons in the cosmic microwave background provide us with a detailed snapshot 
of the conditions of the Universe soon after the Big Bang, photons produced in heavy ion 
collisions represent excellent probes of the "Little Bang" tested at RHIC [1, 2] and LHC 
[3]. The reason is that, like the Universe after recombination, the quark- gluon plasma 
(QGP) created in those experiments is optically transparent, both because of its limited 
spatial extent and because of the smallness of the electromagnetic coupling a EM . Photons 
produced inside of this medium escape with essentially no rescattering, bringing with them 
valuable information about the conditions of the system at the location of their emission 



While it is possible to compute the production rate of such photons in a plasma at small 
a s using ordinary perturbative techniques (see e.g. [6]), the QGP produced at RHIC and 
LHC seemingly behaves as a strongly coupled fluid [7, 8], thus rendering a perturbative ap- 
proach problematic. This non-perturbative nature of the plasma provides a motivation for 
the application of the gauge/gravity correspondence [9-11] 1 to this system. Unfortunately, 
the holographic dual to QCD is not currently known, but the hope is that studying 'rea- 
sonable' (in the sense of [13]) proxies, such as finite temperature N = 4 super Yang-Mills 
(SYM), one might gain insight about the qualitative features of the strong coupling dynam- 
ics of the plasma. Moreover, some quantities, a notable example being the shear viscosity 
to entropy density ratio [14], turn out to be quite model independent and universal. 

1 See [12] for a review of applications to the QGP. 



[4, 5]. 



- 1 - 



The study of photon emission in a strongly coupled N = 4 plasma using holography 
was initiated in [15], where photons coupled to massless quarks in the adjoint represen- 
tation were studied, and later on continued in [16] and [17] where massless and massive 
quarks in the fundamental representation were the electrically charged particles. Photon 
production in a 'soft wall' AdS/QCD model was addressed in [18] and in a charged black 
hole background in [19]. 

The geometry used in those works was spatially isotropic in the gauge theory directions 
and therefore dual to a rotationally invariant plasma. It is now believed that an important 
feature of the real-world QGP is the fact that the pressures along the beam direction and 
along the transverse directions are unequal in the initial stages of the plasma evolution 
and that the system displays a sizable initial anisotropy [20, 21]. A first analysis of the 
electromagnetic signatures in a strongly coupled anisotropic plasma was carried on in [22] , 
where the geometry of [23] was used as dual. An undesirable feature of that model was 
the presence of a naked singularity, that nonetheless allowed the imposition of infalling 
boundary conditions and the definition of retarded correlators. 

In this paper we revisit the computation of photon production in a strongly coupled 
anisotropic plasma using as gravitational dual the type IIB supergravity solution that 
was discovered in [24, 25] building on inspiration from [26]. This solution is static and 
corresponds to a plasma in anisotropic thermal equilibrium, with the force responsible for 
maintaining the equilibrium being a position-dependent theta-angle, as we shall review 
below. Of course, a real-world plasma undergoes a non-trivial time evolution, but the 
assumption of staticity might be motivated by considering time scales much shorter than 
the characteristic evolution scale, which is as long as a few fm/c (see e.g. [27]). 

The geometry of [24, 25] enjoys several nice features that make it an ideal toy model 
for addressing from first principles the effects of a spatial anisotropy on physical observ- 
ables of interest. Such features include regularity of the fields on and outside the horizon, 
asymptotic AdS boundary conditions, and a UV completion with a solid embedding in type 
IIB string theory. We introduce N f <C N c flavors of massless quarks in this background 
and compute the production rate of photons for arbitrary orientations of the photon wave 
vectors and for arbitrary values of the anisotropy. We also compute the electric DC (i.e. 
zero-frequency) conductivity and compare it with the one of an isotropic plasma which 
is kept either at the same temperature or at the same entropy density. We find that the 
conductivity can be larger or smaller than the isotropic conductivity depending on the 
orientation of the photon polarization with respect to the anisotropic direction. 

Other physical properties of the anisotropic geometry considered here have already 
been studied in the literature. These include the shear viscosity to entropy density ratio 
[28, 29], the drag force experienced by a heavy quark [30, 31], the energy lost by a quark 
rotating in the transverse plane [32], the stopping distance of a light probe [33], the jet 
quenching parameter of the medium [31, 34, 35], and the potential between a quark and 
antiquark pair, both static [31, 34, 36] and in a plasma wind [36]. A non-interacting 
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version of the model of [24, 25] has also been considered in [37] and the corresponding 
phase diagram has been obtained. 

The structure of the paper is as follows. In Sec. 2 we derive the field theoretical 
expression for the photon production rate in the presence of anisotropy. In Sec. 3 we review 
the gravitational dual to the anisotropic system we consider and describe the introduction 
of flavor branes in this background. In Sec. 4 we compute the photon production rate and 
conductivity from holography and present our results. We conclude with a discussion and 
an outlook in Sec. 5. 

2 Photon production in an anisotropic plasma 

The gauge theory we shall consider is obtained via an isotropy breaking deformation of 
four-dimensional M = 4 super Yang-Mills (SYM) with gauge group SU(N C ), at large N c 
and large 't Hooft coupling A = gy M N c . The deformation consists in including in the action 
a theta-term which depends linearly on one of the spatial directions, say z, [26] 

Ssu(N c ) = SV=4 + ^2 / TrFAF, 9(z) oc z , (2.1) 

where the proportionality constant in 6{z) has dimensions of energy and will be related to 
the parameter a that we shall introduce in the next section. The rotational 50(3) symmetry 
in the space directions is broken by the new term down to SO (2) in the xy-plane. For this 
reason we shall call the z-direction the longitudinal (or anisotropic) direction, while x and 
y will be the transverse directions. 

This theory has matter fields in the adjoint representation of the gauge group. We can 
also introduce N t flavors of scalars < I> a and fermions ^ a in the fundamental representation, 
with the index a taking values between 1 and N t . With an abuse of language, we will refer 
to these fundamental fields indistinctly as 'quarks'. 

To study photon production we turn on a dynamical photon by including a U(l) kinetic 
term in the action (2.1) and a coupling to the fields that we want to be charged under this 
Abelian symmetry. In order to realize a situation as similar as possible to QCD, we require 
that only the fundamental fields be charged while the adjoint fields are to remain neutral. 
The action for the resulting SU(N C ) x U(l) theory is then 

S = S sum -\jd 4 x (F^ - 4e A» J™) , (2.2) 

where F^ u = d^Ay — d v A^ is the U (1) gauge field strength, e the electromagnetic coupling, 
and the electromagnetic current is given by 

J™ = ^ + V (XV$) - ^ (P^T $ . (2-3) 

with an implicit sum over the flavor indices. The fundamental fields are differentiated using 
the full SU(N C ) x U(l) connection, T>^ = D^ — ieA^, while the 5C/(A ? c )-covariant derivative 
acts on the adjoint matter. 
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Figure 1. Diagrams contributing to the two-point function of electromagnetic currents to leading 
order in the electromagnetic coupling a EM - The external, dotted lines correspond to photons with 
momentum k, while the shaded blobs denote full resummations of the SU(N C ) diagrams, to all 
orders in A. 

We do not know the gravitational dual of the full SU (N c ) x J7(l) theory, but fortunately 
this will not be necessary for our purposes. It was in fact shown in [15] that to compute the 
two-point correlation function of the electromagnetic current (2.3) to leading order in the 
electromagnetic coupling a EM , it is enough to consider the SU(N C ) theory only, whose dual 
is known from [24, 25]. Our computation will then be to leading order in a EM , since the 
coupling of the photons to the surrounding medium is small, but fully non-perturbative 
in the 't Hooft coupling A of the SU(N C ) theory. This is depicted diagrammatically in 
Fig. 1, where the shaded blobs stand for the all-order resummations of the SU(N C ) theory 
diagrams, while the external legs are the U{1) theory photons. 

In general, photon production in differential form is given by the expression [15, 17, 38] 



dT 



dk (2vr) 3 2|A:| 



,$(^ e f s) (fc) £ [ s) W^W 



s=l,2 



k°=\k\ 



(2.4) 



with k^ = (k°, k) the photon null momentum and &(k) the distribution function, which 
for thermal equilibrium as in our case reduces to the Bose-Einstein distribution n#(/c°) = 
l/(e k °l T - 1). The spectral density is x^u(k) = -2 Im G*(k) with 



G* u (k) = -i I d*xe~^ 6(t)<[Jn*), 4 BM (0)]> 



(2.5) 



the retarded correlator of two electro-magnetic currents J* 



Each term of the sum in (2.4) stands for the number of photons emitted with polar- 
ization vector e*( s ). The two polarization vectors are mutually orthogonal and orthogonal 
to k. Given the SO{2) symmetry in the xy-plane, we can choose without loss of generality 
k to lie in the rrz-plane, forming an angle with the z-direction - see Fig. 2. Specifically, 
we set 

k = fc o (sin0,O,cos0). (2-6) 
This means that we can choose the polarization vectors as 



= (0, 1,0) , e (2 ) = (cos 0,0,- sin 0) • 



(2.7) 
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Figure 2. Photon momentum and polarization vectors. Because of the rotational symmetry in 
the xy-jA&ne, the momentum can be chosen to be contained in the xz-plane, forming an angle 9 
with the z-direction. em is oriented along the y-direction and en) is contained in the xz-plane, 
orthogonally to k. 



Production of photons with polarization is then proportional to Xyy ~ ^ m {Jy M Jy M ), 
whereas for those with polarization e^ 2 ) it is proportional to 2 

e (2) e (2) x ^ v = cos2 Xxx + sin2 Xzz ~ 2 cos sin Xxz ■ ( 2 - 8 ) 
From this expression we see that we need to compute the three correlators 

/"fR / tEM tEM\ ^R / tEM tEM\ f~iB. I tEM tEM\ (c, c\\ 

g xx ~ {J x j x ) , g zz ~ \J Z j z } , g xz ~ {J x j z } i^.yj 

and then add them up according to (2.8). In the following section we will see how these 
correlators can be obtained from gravity. 



3 Gravity set-up 

The dual gravitational background for the theory (2.1) at finite temperature is the type 
IIB supergravity geometry found in [24, 25], whose string frame metric reads 

ds z = ^ (-BFdt 2 + dx 2 + dy 2 + Udz 2 + d 4r )+ L 2 e^dn 2 , (3.1) 
u z \ F J 

with T~L = and fis the volume form of a round 5-sphere 

dfi§ = d$ 2 + cos 2 dip 2 + sin 2 dft§ . (3.2) 
2 Note that Xxz = Xxx, see e.g. [15]. 
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Figure 3. Metric functions for a/T ~ 4.4 (left) and a/T ~ 86 (right). 



The gauge theory coordinates are (t, x, y, z) while u is the AdS radial coordinate, with the 
black hole horizon lying at u = u H (where T vanishes) and the boundary at u = 0. As 
mentioned already, we refer to the z-direction as the longitudinal direction and to x and y 
as the transverse directions. L is the common radius of the first metric factor and of the 
sphere and is set to unity in the following. Besides the metric and the dilaton 0, the forms 

F 5 = 4 (J2 6 + *Q 5 ) , F x = adz (3.3) 

are also turned on, with a being a parameter with units of energy that controls the degree 
of anisotropy of the system. The potential for the 1-form is a linear axion x = az - The 
functions B, J 7 , and <j> depend solely on u and are known analytically in limiting regimes of 
low and high temperature and numerically in intermediate regimes [25]. Representatives 
of these functions for two values of a/T are plotted in Fig. 3. For u — > (independently of 
the value of a) they asymptote to the AdS§ x S 5 metric, J- = B = 7~L = 1 and (j) = 0. For 
a = they reduce to the black D3-brane solution 

u 4 

B = n = l, <f> = X = 0, F=l--t, (3.4) 
which has temperature and entropy density given by [39] 

1 7T 2 

T iso = , s iso = — iV c 2 T 3 . (3.5) 

7TU H Z 

The temperature and entropy density of the anisotropic geometry are given by [25] 
m 6-1^^(16 + a 2 ulel^) iV c 2 5, 



16"7TUh 27ru| 



(3.6) 



where 4> H = <j)(u = u H ) and B H = B(u = n H ). As depicted in Fig. 4, the entropy density 
of the system interpolates smoothly between the isotropic scaling above for small a/T and 
the scaling [25, 26] 



i 



s ^3.21iV c 2 T 3 (|) 5 , (3.7) 
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Figure 4. Log-log plot of the entropy density as a function of a/T. The dashed blue line is a 
straight line with slope 1/3. 

for large a/T, the transition between the two behaviors taking place at approximately 
a/T ~ 3.7. 

A feature of the anisotropic geometry of [24, 25] is that some physical quantities (such 
as, for example, the energy density and pressures) depend not only on the dimensionless 
ratio a/T but also on a and T separately. The reason for this is a conformal anomaly that 
appears during the renormalization of the theory introducing a reference scale. Fortunately, 
as we shall see in the following, all the quantities computed in this paper are not affected 
by this anomaly and are only functions of the two independent dimensionless ratios that 
can be built out of a, T and ko (e.g. a/T and ko/T). 

3.1 Flavor branes 

The introduction of N f flavors of massless quarks is achieved by placing N t probe D7-branes 
in the background (3.1). Given that the quarks are massless, any temperature different 
from zero in the background will imply that these quarks are in a deconfined phase (the 
'black- hole embedding' [40]) and also that the induced metric over the D7's will be the 
metric (3.1) with a three-sphere instead of a five-sphere. The complete system can then 
be thought of as a D3/D7 system with two different kinds of D7- branes, one kind sourcing 
the anisotropy [25, 26] 3 and the other kind sourcing flavor [41, 42]: 
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(3.8) 



These branes are smeared homogeneously along the ^-direction and can be thought of as giving rise to 
a density nm = N-or/Lz of extended charges. This charge density is related to the anisotropy parameter a 
through a = gY M nD7/47r [25]. 
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Throughout our calculations, the value of N m (and consequently of a, see footnote 3) can 
be arbitrary, but we will assume that flavor is quenched, N f <C N c . As explained in [25], the 
full backreaction of the 'anisotropic' D7-branes is already included in the geometry (3.1), 
through the presence of the axion field x- I n the following, when talking about D7-branes 
we will refer exclusively to the flavor D7-branes, which are probes in the fixed background 



As argued in [15], to study photon emission to leading order in e it suffices to evaluate 
the correlators needed for (2.4) in the SU(N C ) gauge theory with no dynamical photons. At 
strong 't Hooft coupling and large N c , these correlators can be calculated holographically. 
Indeed, global symmetries of the gauge theory are in one-to-one correspondence with gauge 
symmetries on the gravity side, and each conserved current of the gauge theory is dual to 
a gauge field on the gravity side. 

Let A m (m = 0, . . . , 7) be the gauge field associated to the overall U (1) C U(N { ) gauge 
symmetry on the D7-branes. Upon dimensional reduction on the 3-sphere wrapped by 
the D7-branes, A m gives rise to a massless gauge field (A^^Ay), three massless scalars, 
and a tower of massive Kaluza-Klein (KK) modes. All these fields propagate on the five 
non-compact dimensions of the D7-branes. We will work in the gauge A u = 0, 4 and we 
will consistently set to zero the scalars and the higher KK modes, which are not of interest 
here. The gauge field A^ is the desired dual to the conserved electromagnetic current J® M 
of the gauge theory. According to the prescription of [10, 11], correlation functions of J^ M 
can be calculated by varying the string partition function with respect to the value of A^ 
at the boundary of the spacetime (3.1). 

We proceed now to writing down the action for the D7-branes. It is easy to realize 
that there is no Wess-Zumino coupling of the branes to the background F5, because of the 
particular brane orientation that has been chosen. A priori there could be a coupling to 
the background axion 



but this would be quartic in the U{1) field strength F = dA and therefore irrelevant for 
the present computation, since only 2-point functions are needed. 

This means that the Dirac-Born-Infeld (DBI) action is all we need to consider: 



where g is the induced metric on the D7-branes and T D7 = 1/ l {2ir£ s )' ' g s l s is the D7-brane 
tension. To obtain the equations of motion for A^ it suffices to expand the action above 
and use the quadratic part only 



4 This gauge choice will be immaterial in the following, since we shall switch to gauge invariant quantities, 
but it has the advantage of simplifying our formulas. 



(3.1). 




(3.9) 




(3.10) 




(3.11) 
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where F 2 = F mn F mn . The embedding of the branes inside the of the geometry can be 
parametrized by an angle cos i? = ip(u) an d the induced metric on the branes is given by 

dsl 7 = -A-TB dt 2 + dx 2 + dy 2 + U dz 2 ) + l ~^ltf Fe *y 2 du 2 



1 



+e2*(l-V 2 )^3. (3.12) 

The case of massless quarks corresponds to the equatorial embedding of the D7-branes, 
namely to ip = 0. After the dimensional reduction on the three-sphere, the action reduces 
to 

S = -K D7 / dtdxdu ^F 2 , (3.13) 



where 

K-D7 = 2ir A N c T m £i = ^ N ° N t > ( 3 - 14 ) 
and F m is restricted to the components m = (yu, u). 

As argued in [17, 40], in order to calculate the photon emission rate, we may consis- 
tently proceed by assuming the equatorial embedding of the D7-branes to be fixed in the 
absence of the gauge field, and then solving for the gauge field on that embedding. No 
modes of the metric or of the background fields will be excited. As mentioned above, we set 
to zero the components of the gauge field on the three-sphere wrapped by the D7-branes 
and we Fourier decompose the remaining components as 

^^e- ik ° t+ik - x A^k ,k,u), fc = fc o (sin0,O,cos0). (3.15) 

This is possible because the state we consider, albeit anisotropic, is translationally invariant 
along the gauge theory directions [25]. Doing so, the equations for the gauge field deriving 
from (3.13) split into a single, decoupled equation for A y , which reads explicitly (primes 
denote derivatives with respect to u) 

(Mg uu g yyA' y )' - Mgyy {g tt k 2 + fk 2 x + g zz k 2 z ) A y = 0, (3.16) 

plus a coupled system of three equations for the remaining components At >x ,z' 

(Mg uu g tt A' t y - Mg tl [g xx k x (k x A t - k A x ) + g zz k z (k z A t - k A z )} = , (3.17) 
(Mg uu g xx A' x Y - Mg xx [g tt k (k A x - k x A t ) + g zz k z (k z A x - k x A z )] = , (3.18) 
(Mg uu g zz A' z )' - Mg zz [g tt k (k A z - k z A t ) + g xx k x (k x A z - k z A x )] = . (3.19) 

The inverse metric can be read off directly from (3.1) and M is the factor appearing in the 
action (3.13): 

M ee —f- . (3.20) 

These are the set of equations that need to be solved (with the appropriate boundary 
conditions), in order to obtain the correlation functions of the electromagnetic currents 
jem_ ^y e ^ ^ n following section. 
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4 Photon production from holography 



The Lorentzian AdS/CFT prescription to compute the correlation functions of the elec- 
tromagnetic current is given in [43] (see also e.g. [44-46]). We need to isolate the terms 
with two radial derivatives in the action (3.13), write the latter as a boundary term by 
integration by parts and compute it on-shell. Varying with respect to the values of at 
u = (which correspond to the boundary theory currents) gives the desired correlation 
functions. 



We start by writing the boundary action as 



S e = -2K D7 J dt dx e (-J^ApAl + A X A' X + A y A' y + e^A z A')j 



, (4-1) 



where the limit e — > is intended. 



From this boundary action and the equations of motion (3.16)-(3.19) we see that we 
can carry on the calculation for the spectral density Xyy independently of all the others 
and the calculation is very similar to the one in [17]. The remaining spectral densities in 
(2.8) are more involved and will require a strategy that was originally developed in [47] to 
deal with operator mixing. 



4.1 Isotropic limit 

We shall see in the following that in order to properly compare our results with the photon 
production rate in an isotropic plasma we will need expressions for the isotropic correlators 
both as a function of the temperature of the plasma and as a function of the entropy density 
of the plasma. 

We start then by considering the isotropic limit, which is obtained by plugging (3.4) 
in the equations above. We may also set 9 = tt/2, since in the isotropic system we have an 
50(3) symmetry that permits to align the photon momentum along a particular direction, 
say the x-direction in our choice. The equation to solve is obtained from (3.16) and reads 5 

- » ( l - ff< *3 AL + kl i A - ■ (4 - 2) 

The solution of this equation that is infalling at the horizon reads 
t» to 

( « 2 V' T ( v?\~^ I 1 + i 1+i "l-u 2 \ , . 

A ~ = I 1 - <J ( 1 + ^J ^ (' - -1— ■ —j—- 1 - it0; isr ) ' (4 ' 3) 



where the dimensionless frequency 



ko k Q u H 

to = — — = (4.4) 

2vrT 2 v 1 



5 This equation can be seen to be precisely equal to eq. (4.18) of [17] by doing the change of coordinates 

Where = «H \Z2u t here ■ 
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has been introduced, as customarily in the quasi-normal mode literature. 



Rotational invariance implies that the sum over polarizations in (2.4) be equal to the 
following expression 



Xiso = -21m (Gy y + G* ) ;so = -41m Gf ao , 
which can be computed from the prescription in [43] using 

4K D7 



lim ( 1 



|A so (t», 0)| 2 u^o 



U7, I U 



-A* so (™,u)A[ so {w,u). 



(4.5) 



(4.6) 



The imaginary part of this expression turns out to be independent of u and it is actually 
easier to evaluate the limit at the horizon rather than at the boundary. The result of the 
limit is 

«)=i^. (4.7) 



lim ( 1 

u-)-mh 



u 4 \ 1 



Putting everything together, one finally obtains that the spectral density is given by [17] 

8Af UT xx> 



Xiso 



2"> \ 2 Fi (1 - i±i», -ifiw, 1 - in; \) \ 



where 



Mr, 



K m = 1 

: ul 4 



N r N r T 2 . 



(4.8) 



(4.9) 



The dimensionless combination Xiso/8A/" D 7tx> is plotted in Fig. 5. We can now express Ad 




5 10 15 

Figure 5. The isotropic correlator Xh 



20 w 



and to in the formula above either in terms of the temperature or in terms of the entropy 
density of the isotropic plasma, using (3.5). We obtain 



N c N t Tk Q 



227TT7T 



(4.10) 
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and 



Xiso(s) = - i 5 — — —2 , Ws = 4 i 1 • I 4 - 11 ) 

2 TOs -37r3 | 2 Fi (1 - ^rx> s , -^»„ 1 - it» s ; \) | 23vr3 S 3 

In the following sections we shall use these expressions to normalize the anisotropic results. 
The two different normalizations will allow us to compare the anisotropic photon production 
either with an isotropic plasma at the same temperature but different entropy density or 
with an isotropic plasma at the same entropy density but different temperature. 



4.2 Spectral density for the polarization en\ 

After the isotropic case warm-up, we proceed to compute Xyy This can be obtained 
independently of all the other spectral densities, given that (3.16) is not coupled to the 
remaining equations of motion and that there is no mixing term in the action between A y 
and the other gauge fields. The corresponding correlation function is 

4K 

G vv = i a fu nM2 ^Q(u)A;(k ,u)A> y (k 0l u) , (4.12) 
\A y (k ,0)\ z 



where 



The spectral density reads 



Qiu) = — . (4.13) 



N Nf , 
X(l) = Xyy = 2 / 2 Im l[m n Q( u ) A l(. k o, u)A y (k , u) . (4.14) 

In the anisotropic case the metric components appearing in the equations of motion 
are only known numerically for generic values of a/T. This means that we have to resort 
to numerics in order to integrate (3.16) and obtain the solution for A y . To impose the 
appropriate boundary conditions at the horizon we use the near-horizon expansion of the 
metric given in [25] and expand (3.16). We see that close to the horizon A y behaves like 
(u — u s ) u a y (u), where a y (u) is some function regular at u H . The exponent v turns out to 
be given by the usual result 

* = ±V »=^, (4-15) 

where the temperature appearing in this formula is now the anisotropic one, given by 
(3.6). To impose the infalling wave condition at the horizon and therefore get the retarded 
correlator, we choose the negative sign for v. Once the infalling solution has been selected, 
the expansion is unique up to an overall multiplicative factor. 

Now we carry out the numerical integration of (3.16) for a number of values of the 
anisotropy and of the angle 9 formed by the momentum k and the z-axis. The spectral 
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Figure 6. Plots of the spectral density X(i) corresponding to the polarization em normalized 
with respect to the isotropic result at fixed temperature Xiso(T). The curves correspond from 
top to bottom to the angles 6 = 0, 7r/8, 7r/4, 3tt/8, tt/2. The four plots correspond to the cases 
a/T = 1.38 (a), 4.41 (b), 12.2 (c), 86 (d). 

density Xyy is plotted in Figs. 6 and 7, where it is compared with an isotropic plasma at 
the same temperature or at the same entropy density. Note that Xyy i n the isotropic limit 
goes to half of x iso , since the latter contains the (equal) contribution from Xzz as well, see 
(4.5). 

Note also that at large ko (more specifically, for ko 3> a) the spectral density does not 
converge to the isotropic limit (the horizontal line at 1/2). This is perhaps surprising at 
first sight, but it is explained remembering that our results depend on two independent 
dimensionless ratios. Sending k$ — > 00 while keeping a and T fixed is not the same as 
sending a — > with ko and T fixed, since the ratio a/T remains finite in the first limit. As 
a consequence, the large-feo limit of our result will reproduce the isotropic limit only if we 
send a/T — > at the same time. 6 



We thank David Mateos for comments on this point. 
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Figure 7. Plots of the spectral density X(i) corresponding to the polarization em normalized 
with respect to the isotropic result at fixed entropy density Xiso(s). The curves correspond from 
top to bottom to the angles 6 = 0, 7r/8, 7t/4, 37r/4, tt/2. The four plots correspond to the cases 
aAr 2/3/ s i/3 = o 80 ^ 2.47 (b), 6.24 (c), 35.5 (d). 

The zero-frequency limit of the spectral density gives the electric DC conductivity. 
For photons with polarization e(i) this would be the conductivity along the transverse 
y-direction. The quantities 

a a) (T) = lim X(1 > = lim 2^%- , 

a (l) (s) = lim X(l \ . = lim (4.16) 

k ' k ^0 X(l),iso(s) *o-+0 Xiso(s) 

as functions of a/T and of aN^/ 3 / s 1 / 3 are reported in Fig. 8. The plot on the left of 
that figure coincides with the perpendicular conductivity a± obtained in [28]. Note that 
the conductivity is to be normalized with respect to the isotropic result with the same 
polarization, i.e. lim a _>o = X(i),i so m the present case. 
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4.3 Spectral density for the polarization e/ 2 ) 

Having found the correlator corresponding to the photon polarization em , we proceed now 
with the correlators corresponding to e^) ■ 

The procedure to find them is clearer if done in terms of the gauge invariant fields 
Ei = 8iAq — doAi. Equations (3.17)-(3.19) can be rewritten in terms of Ei with the aid of 
the constraint 

- ~A' t + smOA' x + cos6e*A' z = 0, (4.17) 
resulting in the set of equations 



E" + 



E" + 



log 



log 



e 4-. 



BT 



u 



+ |(io g ^y 



K + yE x + 



?^k x k z 
k 2 



(log&F)% = 0, 



(4.18) 



e4 ■ 



BT 



+ ^£(iog e ^y 



E> z + k yE z + k -^{\o g e*BF)>E> x = V. 



(4.19) 



where k 2 = — k 2 — e^k 2 . Notice that for 6 = (momentum along the z-direction) and 
9 = 7r/2 (momentum along the x-direction) the equations for E x and E z decouple. 



The action (4.1) at u = e can be written in terms of these fields as 



-2K r 



I 



dt dx 



e'^y/BT f fk. 



U fcn k 2 



KB ^-e*k 2 )E x E x + k 2 E y E' y 



+e^k x k z {E x E z )' + 



( 



v - kl ) e*E z E> z 



(4.20) 
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To perform the variations of the action, it is important to notice that differentiating with 
respect to Ai =X) y tZ is just ko times the result of differentiating with respect to Ei, and to 
differentiate with respect to Aq we have to sum ^ ki g^r . 

Since we need to take the limit e — > let us compute the on-shell value of the integrand 
in (4.20) using the near-boundary expansion of the metric. This expansion can be found 
in [25] and is given by 

T= l + ^aV+V + ^aVlog« + 0( M 6 ), (4.21) 

B = 1 - ^a 2 « 2 + S 4 u 4 - ^a 4 u 4 logu + 0(« 6 ), (4.22) 

a 2 2 / 1152234 + 121a 4 \ 4 a 4 , fiN . . 

4> = --u 2 + ( ^_ J u 4 - -u 4 log u + 0(u 6 ), (4.23) 

where T± and B\ are coefficients undetermined by the asymptotic equations of motion that 
can nonetheless be extracted from the numerical solution. They turn out to be irrelevant 
for our purposes. 

Using (4.21)-(4.23) we can solve equations (4.18) and (4.19) perturbatively to find 
E, - *f > + EP cos*.* - 1 (^ co^ + Hff>) + O (» 6 ) , (4,4, 



E z = - 4 2 ) sin^ 2 + £i 4 V - C ° s2 i g 6 a2fc ° 2 (4°) tanfl + « 4 logu + O 



(4.25) 

where E$\ E?\ E%\ and 2# } are undetermined coefficients of the expansion. They can 
be extracted from the numerical solution. Using the former expressions we can write the 
value of (4.20) close to the boundary as 

S e = -2K Dr J dtdx [d + C 2 + C 3 + . . . + O (u 2 )] u=t , (4.26) 

where 

d = -- sin 2 9 4°) 2 - l - cos 2 9 £(°) 2 - cos 9 sin 9 E^E® , 
1 



o 



1 + omsSfl^pj + 48 ^ ^(0)^(4) _ 10sin ^(0)^(2) + 



cos6 ~ x ' a 2 x ~* ~ z ~ x ' a 2 ' 

C a = - ^i 0) sin 9 + £^ 0) cos log u , (4.27) 

and the ellipsis stands for the terms in the y-components that have already been dealt with 
in the previous section. 

We notice that £3 diverges as we take the u — > limit, which implies that the boun- 
dary action is sensitive to the anomaly of the background found in [25]. In principle we 
would need to holographically renormalize this action, but the key observation here is that 



- 16 - 



according to (2.8) the contribution of this divergent term to the production of photons 
with polarization e*( 2 ) is proportional to 

5 2 C b 2 C b 2 L 
cos 2 8 ^+sin 2 ^-2cos0sin0 — m 3 m = 0, (4.28) 

6E^ 2 6E^ 2 SE^SE^ 

and therefore vanishes identically. This means that we only need to consider the contribu- 
tions to the correlators coming from C\ and £2, which are finite and turn out to be given 
by 



Gxx 3 ,2 ■ 2a , l + 5cos20<^ 2) , 16tan0^i 4) 

- -k 2 S m 2 9 + Q _ Q ,_. + (4.29) 



-4K D7 ~ 4 3cos0 SE <S>) a 2 6E jfl) ' 



G% z l,o , „ lOsin0^i 2) . IQSE^ 



-k 2 cos 2 8 - ^ + ~ ^ , (4.30) 



-4^dt 4 u 3 ^(0) a2 §E I 

G% ;2 . . . , l + 5cos20<^ 2) 

— = - ki cos 9 sin 9 H 7^ 

-2K m 3cos6» 5E f) 



16 tan 6> <5£| 4) 10 sin 8 5E X 2) 16 5E {4) 



+ 70— k- ( 4 - 31 ) 



« 2 <^ 0) 3 5jB (0) «2^(0) 

Summing up these expressions according to (2.8) one finds the simple result 

X(2) = ^2) e (2) = 16i^D7lm ( cos °— Joj ~ sind J^{o)J ■ ( 432 ) 

To determine how the coefficients E^ and E^ vary with respect of E x °^ and E^ it 
is convenient to write the solutions E x and E z as a vector. When solving the equation of 
motion for 

E-(*) (4.33) 

near the horizon by a power series expansion (u — u^Y we get that the exponent v for both 
components of this vector is the same as that for the A y mode, namely v = ±ivo/2. After 
imposing the infalling wave condition, the rest of the power series is linearly determined 
by the value of E at the horizon. Integrating from the horizon using any choice of such a 
vector would pick a particular solution to the equations of motion. These are linear and 
therefore their general solution can be written as a linear combination of any two linearly 
independent solutions Ei = (E x> i, E Zt \) T and E2 = (E x> 2, E Z ^) T ■ 

Since what we need to know is how the solution depends on E x ^ and E^ , the relevant 
base is the one formed by the two particular solutions that at the boundary reach the values 

(4.34) 





because in terms of this base the general solution can be written as 

E = Ei + Ef ] E 2 . (4.35) 
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Close to the boundary the general solution E will be described by the expansions (4.24) 
and (4.25). Prom (4.35) we see that its corresponding coefficients will be given by 

4 2) = 4 0) E% + E?> E% , E?) = 4 0) E<$ + Ef) E% (4.36) 

where Ef] and Ef] (with i = 1,2) are the coefficients for the expansions of Ei and E2 
close to the boundary. Prom (4.36) it follows, for instance, that 



wi 2) _ m SEP _ m 

~m ~ E ^ 7pm ~ E ^ ' (4 - 37) 



8E£> ' 5E\ 

which can now be used in (4.32). In practice, since we can only solve the equations of motion 
numerically, we find any two linearly independent solutions, say E a = (E XiCl , E z ^ a ) T and 
Eb = (E xb , E z b ) T , to construct Ei and E2 as the linear combinations 

r(°) E -£ (0) Ek E {0) Ek - E {0) E 

z b a ^z.a -C^b J^x,a *-> x b - c -'a 

El = E m E {0) - E {0) E {0) ' ^ = E {0) E (0) - E {0) ' 

whence we can see that 

P (0) p,(2) p(0) „(2) „(o) p (2) „(o) e,(2) 

F (2) _ ^z,b ~ & x fi ( 2 ) _ ~ ^x,b „ q s 

^ ~ p(0) p(0) F (0) „(0) ' ^.2 - p (0) „(0) F (0) p (0) • 

With this groundwork in place, we are now in the position of using these expressions to 
numerically obtain the coefficients we need to input in (4.32) to find X(2)- The results are 
plotted in Figs. 9 and 10, where we compare, for different angles 9, the photon production 
with polarization e/ 2 ) with the one of an isotropic plasma at the same temperature or at 
the same entropy density. The 9 = curve for \(2) in ai l the plots matches the 9 = curve 
for X(i)- This is expected, since in this case equations (3.16) and (3.18) are identical, being 
gxx _ g-yy^ Q n th er hand, one does not expect a similar match between the 9 = tt/2 
curve for %(2) an d the 9 = tt/2 curve for X(l)i since in this case equations (3.16) and (3.19) 
are different, being g yv ^ g zz , and indeed we get different curves. The comment we had 
f° r X(i) about the fact that the large-fco spectral density will reproduce the isotropic limit 
only if a/T — > also applies here. 

For photons with polarization along e( 2 ) the conductivities 

a m {T) = lim X ® = Hm 2 ^%- , 
( )V ' *o->0 X(2),i SO (T) fc ->0 Xiso(T)' 

a (2) (s) = lim X{2 \ . = lim 2^T ( 4 - 4 0) 

1 ' fco^O X(2),iso(s) ^0^0 Xiso(s) 

depend not only on the anisotropy a, as was the case for the polarization along the y- 
direction, but also on the angle 9. Plots of ci(2) as a function of a/T and aN^ 3 /s 1 / 3 (for 
fixed 9) and as a function of 9 (for fixed a/T) are reported in Figs. 11 and 12. For 9 = 
the conductivity is along the transverse rr-direction (i.e. the direction of the polarization 
vector) and we recover the same results as (4.16) and a±_ of [28]. Similarly, for 9 = tt/2 the 
conductivity is along the longitudinal z-direction and we recover a z of [28]. 



- 18 - 




Figure 9. Plots of the spectral density %(2) corresponding to the polarization £( 2 ) normalized 
with respect to the isotropic result at fixed temperature Xiso(T). The curves correspond from 
top to bottom to the angles 6 = 0, 7r/8, 7r/4, 3tt/8, tt/2. The four plots correspond to the cases 
a/T = 1.38 (a), 4.41 (b), 12.2 (c), 86 (d). 

4.4 Total production rate 

We have now all the ingredients to plot the total production rate (2.4). We convert this 
quantity to emission rate per unit photon energy in an infinitesimal angle around 9. Using 
that the photon momentum is light-like we get 

2a EM N c N t T 3 d(cos9) dk° ~ 2N c N t T 2 e 2 ^> - 1 + X{2) > ' ( ' ' 

This quantity for different values of a/T and for different values of 9 is plotted in Fig. 13. 
The isotropic result at the same temperature is [17] 

-1 dT lan w 2 



2a EM N c N t T 3 d(cos9) dk° 2 ™ (e 2 ^> - 1) \ 2 F t (l - l±*w a , -^to s , 1 - *w a , 



_!+»«« 1 _ ^ ■ IM 2 ' 

(4.42) 

and is plotted as a black dotted curve in the figures. 
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Figure 10. Plots of the spectral density X(2) corresponding to the polarization e/2) normalized 
with respect to the isotropic result at fixed entropy density Xiso(s)- The curves correspond from 
top to bottom to the angles 6 = 0, 7r/8, 7r/4, 37r/4, tt/2. The four plots correspond to the cases 
aA r2/ 3 / s i/3 = o 80 ( a ) ; 2.47 (b), 6.24 (c), 35.5 (d). 



5 Discussion 



In this paper we have computed the photon production rate and the electric conductivity 
of a M = 4 anisotropic plasma modeled holographically using the geometry of [24, 25]. The 
plasma is infinitely extended, in thermal equilibrium, strongly coupled, has a large number 
of colors N c and a quenched number of flavors, N { -C iV c . It is rotationally invariant in 
the x- and y-directions, but not in the z-direction, which we have called the anisotropic 
direction. In a heavy ion collision experiment, the latter would correspond to the beam 
direction, whereas the xy-plane would correspond to the transverse plane. We have con- 
sidered arbitrary orientations of the photon momentum with respect to the anisotropic 
direction, as well as arbitrary values of the anisotropy parameter. 

The spectral densities are sensitive to the degree of anisotropy of the system and to the 
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Figure 11. Plot of the conductivity CT(2) corresponding to the polarization ep) as a func- 
tion of a/T (left) and of aN^ 3 /s 1 / 3 (right). The curves correspond from top to bottom to 
(9 = 0, tt/8, tt/4, 3tt/8, tt/2. 
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Figure 12. Plot of the conductivity ct( 2 ) corresponding to the polarization €( 2 ) as a function of 0. 
We plot the cases a/T = 1.38 (black, solid curve), 4.41 (red, dashed curve), 12.2 (magenta, dotted 
curve), 40 (green, very coarsely dashed curve), 86 (orange, dot-dashed curve), 120 (brown, coarsely 
dashed curve). 



angle between the anisotropic direction and the photon wave vector, with the sensitivity 
increasing with the photon energy. This increasing sensitivity has also been found in 
anisotropic plasmas at weak coupling, see e.g. [48]. The vo — > limit of the spectral densities 
(i.e. the DC conductivity) is independent of the angle 9 for photons with polarization em 
(along the transverse y-direction) and dependent on 9 for photons with polarization e^) 
(in the :rz-plane), as expected on general grounds. 

A generic feature emerging from our study is that the spectral densities (and conse- 
quently the total production rate) are higher for photons with longitudinal wave vectors, 
9 = 0, than for the ones with transverse wave vectors, 9 = tt/2. The enhancement along 
the longitudinal direction is perhaps not surprising, since that is the preferred direction in 
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Figure 13. Plots of the total production rate. The curves correspond from top to bottom to the 
angles 9 = 0, tt/8, it/A, 37r/8, tt/2. The four plots correspond to the cases a/T = 1.38 (a), 4.41 (b), 
12.2 (c), 86 (d). The temperatures in the four cases are, respectively, T = 0.32, 0.33, 0.36, 0.49. The 
isotropic result at the same temperature is the dotted black curve. In (a) the curves are virtually 
indistinguishable from the isotropic one. 

the expansion of the plasma. This enhancement also takes place in the electric conductiv- 
ity, which is larger (smaller) than the isotropic value along the perpendicular (transverse) 
direction, irrespectively of whether the comparison is made with an isotropic plasma at 
the same temperature or at the same entropy density. Note that the conductivity is in the 
direction of the electric field, which coincides with the direction of the polarization vector 
and not of the photon momentum. 

The source of anisotropy in [24, 25] is a position-dependent theta-angle in the gauge 
theory, or, equivalently, a position-dependent axion in the gravity background. One might 
wonder how universal are the behaviors we have obtained, i.e. how they depend on the 
source of anisotropy used. We notice that the equations of motion for the gauge field on the 
branes' world-volume are solely determined by the metric and the dilaton. This implies 
that any other source of anisotropy that gives rise to a qualitatively similar metric and 
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dilaton and no -B-field will produce qualitatively similar results for the photon production 
rate, irrespectively of the other background fields (in particular, irrespectively of whether 
or not there is a background axion turned on). 

Photon production in the strongly coupled anisotropic plasma of [23] has been worked 
out in [22]. The geometry in those papers presents a naked singularity, which nonetheless 
allows for infalling boundary conditions and the definition of a retarded Green's function 
in the dual gauge theory. The small frequency behavior of the spectral densities found in 
[22] is totally different from the one found here, due perhaps to the absence of a bona fide 
horizon. In particular, they find that for any non-vanishing anisotropy, even if very small, 
the spectral densities go to zero faster than vo, resulting in a vanishing DC conductivity. 
At large frequencies, on the other hand, their anisotropic results smoothly approach the 
isotropic limit when the anisotropic parameter is sent to zero, as in our case. They find 
that the photon production rate is enhanced (suppressed) with respect to the isotropic 
one in the forward direction and suppressed (enhanced) in the transverse direction if the 
anisotropy is prolate (oblate). This is to be contrasted with our results, where the rate is 
always enhanced with respect to the isotropic one, although it is generally larger in the 
forward direction than in the transverse one. 

The geometry considered in the present paper is static and the corresponding plasma in 
thermal equilibrium. This fact renders a comparison with a real- world plasma problematic, 
at least for small vo, which corresponds to large times. The plasma evolves in fact from an 
initial out-of-equilibrium state toward a final isotropic state and the static approximation 
necessarily breaks down. The inclusion of time evolution in the present context would 
require a time-dependent geometry and a distribution function more general than the 
Bose-Einstein one. Finding such solution might be somewhat simplified using the ideas of 
[49], but it is beyond the scope of the present paper. 

It would be interesting to extend our analysis by including a non-vanishing mass for 
the quarks, thus bringing the model closer to QCD. A non-zero quark mass is achieved by 
considering embeddings of D7-branes that, unlike the case studied here, are not equatorial. 
The asymptotic distance between the D7-branes and the equatorial plane would determine 
the mass of the quarks in the dual gauge theory. 7 In the isotropic case the photon pro- 
duction rate showed a dependence on this mass [17], so it is clear that such a dependence 
should also be expected in the anisotropic case. In particular, as the mass becomes large 
compared with the temperature or anisotropy, it should be possible to identify the appear- 
ance of highly localized resonances in the spectral functions, indicating the reconstitution 
of the mesons that were melted in the plasma [17]. Understanding the role played by the 
anisotropy in this transition is certainly something worth exploring. 

Finally, another possible extension would be the computation of the production rate 

7 In a strongly coupled plasma there are no quasi-particles and the mass mentioned here would in fact 
simply be a microscopic parameter on which the physics depends, but without the direct interpretation of 
mass of a quasi-particle. 
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of dileptons, which are also known to be a good probe of the properties of the emitting 
medium. 
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